arXiv:1501.07572vl [physics.gen-ph] 1 Dec 2013 


Stability of 2-body orbits 
in retarded gravitation theory (RGT) 

C. K. Raju 

AlBukhary International University 
05200 Alor Setar, Malaysia 
ckr@ckraju.net, ckraju@aiu.edu.my 

January 30, 2015 


Abstract 

The recently formulated retarded gravitation theory (RGT) ex¬ 
plains the non-Newtonian velocities of stars in spiral galaxies, without 
any new hypothesis, and may hence be tested even in the laboratory. 
However, doubts have been expressed that those higher rotation veloc¬ 
ities in RGT may be due to instabilities. We resolve these doubts by 
solving the full functional differential equations of RGT for a model 
2-body planetary system. The solution is stable and closely agrees 
with the Newtonian solution for this planetary case. Thus, the big 
difference between RGT and Newtonian gravity for a spiral galaxy is 
not due to any instability in RGT. 

Keywords: Retarded gravity, galactic rotation curves, functional differential 
equations, experimental tests of gravity 
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1 Introduction 

1.1 Preliminaries 

Recent evidence (from supernovae PJ |2], and the cosmic microwave back¬ 
ground PHI ElE]) indicates that the expansion of the cosmos is accelerating. 
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This has belied long-standing expectations, based on Newtonian gravity or 
general relativity theory (GRT), that the expansion must decelerate. This 
acceleration can, of course, be explained by populating the cosmos with a 
mysterious invisible substance which exerts negative pressure—dark energy. 
However, this has also led to doubts about the validity of our current under¬ 
standing of gravitation. 

Another departure from Newtonian expectations has long been known: 
the rotation velocities of stars in spiral galaxies increase as we move out 
to the edge of the galaxy, instead of decreasing, as expected on Newtonian 
gravity. [7j This departure too has been explained by populating the galaxy 
with invisible dark matter. A further fact contrary to Newtonian expecta¬ 
tions is that the observed rotation velocities of stars reach a constant value 
at the edge of the galaxy. [8] But this too can be explained by advancing 
yet another hypothesis that the dark matter is distributed in the form of a 
halo, with its density peaking where the luminous matter thins out to noth¬ 
ing. Furthermore, this is not just some unseen mass, but specihcally cold 
dark matter—involving some hypothetical form of particles for which there 
is no laboratory evidence. This has led to the skepticism that the “standard 
model” of cosmology “ACDM... relies heavily on two potentially fictitious 
invisible entities”. [9] A theory ceases to be scientific (refutable) if it counters 
every observational surprise with a new hypothesis, and its criticism is inhib¬ 
ited on grounds of “standardised” community opinion. The other common 
alternative is modihed Newtonian dynamics (MONO) [TUI 1^ which modihes 
the gravitational force at scales beyond a “characteristic scale”, but in a 
purely phenomenological way, and goes on to accumulate further hypothe¬ 
ses. 

This suggests the need for a theoretically rigorous approach, where such 
speculative hypotheses (whether about dark matter or about scale-driven 
changes in the gravitational force) are avoided. Accordingly, a new explana¬ 
tion has recently been advanced, m for galactic rotation curves, without any 
new hypothesis. The non-Newtonian velocities in spiral galaxies may be just 
a consequence of Lorentz covariance (which itself is a theoretically essential 
correction [12] to Newtonian physics). 

This seems astonishing at first sight, for special relativistic (-) effects are 
believed to be relevant only at relativistic velocities u ~ c. However, while 
that is dehnitely true for the one body problem, does it remain true for a 
billion-body problem? Even when u -C c could a systematic ^ effect become 
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significant when snmmed over a billion bodies or more, snch as the stars in 
a spiral galaxy? 

This qnestion cannot be answered nsing GRT (or similar geometrical 
theories), since with GRT one is barely able to do the 2-body problem after 
a century. Nor can it be answered with the parametrized post-Newtonian 
(PPN) formulation (whether regarded as an approximation to GRT or as an 
independent theory). For a billion-body problem, Newtonian gravitation was, 
until recently, the only available option, and that obviously cannot be used 
to answer the question whether - effects can accumulate when the number of 
bodies is large. However, the question can now be answered (and the answer 
is yes) within retarded gravitation theory (RGT).[TT] 


1.2 Retarded gravitation theory 

RGT is a Lorentz covariant theory of gravitation which agrees with New¬ 
tonian gravitation in the static case. It was recently formulated as follows. 
Gonsider two particles, with world lines given by Ti(si) and ¥ 2 ( 32 ), where 
Yi and ¥2 are 4-vectors, and si and S 2 are the respective proper times. The 
equations of motion in RGT are 


d^y\ „ 


d‘^y2 ^ 


( 1 ) 


where mi and m 2 are the respective rest masses of the two particles, and 
Fi 2 the 4-force exerted by particle 2 on particle 1 is given by the Lorentz 
covariant expression 


Fi2 — 


kc^ 


(-^2 ret • R 2 ret 

fee" p 

iR2.v2r^" 


_p I_ fec^ (-R2ret-hl) .^ 

)3 {R2ret.V2ret)HV2ret.Vl) 


kc^ (R 2 .R 1 ) 
iR2.V2riV2.Vl) ' 


2 ret 


( 2 ) 

(3) 


Here, k = Gmim 2 , G is the Newtonian gravitational constant, c is the speed 
of light, R 2 ret = h 2 ret “ hi IS the retardation vector, Vi = ^ and V 2 = 
^ denote the respective 4-velocities, and, in ([^, [] 2 ret indicates that the 
quantities with subscript 2 are to be evaluated at the corresponding retarded 
proper time, as explicitly indicated in ([^. The other force F 21 is given by 
interchanging 1 and 2 in ([^. 
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In coordinates, if Yi = (ct,yi(t)), and I 2 = {ct,y 2 (t)), the retarded 
coordinate time ti 2 , in the force Fu acting on Yi at time to, is the root of 
the equation 

c\t-tof=rl,^{y2{t)-y,{to)f, (4) 

satisfying t < to- That is, it is the value of t at the spacetime point where 
the backward null cone from yi(to) intersects the world line Y 2 . The corre¬ 
sponding distance ri 2 is the retarded distance from particle 1 to particle 2. 
A similar equation holds for ^ 21 , the retarded coordinate time in F 21 , the 
asymmetry being only in the arguments of yi and y 2 - 

For non-relativistic velocities, the force in RGT approximately reduces to 
the simpler expression 


T"i2 
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I R2ret 1^2 ret A 

^^^12 c y ’ 


(5) 


which still differs from the Newtonian gravitational force by its (tiny) velocity 
dependence. (It is exactly the Newtonian gravitational force, if the two 
masses are relatively at rest.) 

In a spiral galaxy, all stars in the neighbourhood of a given star are 
rotating in the same direction. Though the rotational velocities are non- 
relativistic (); ~ 10“^ 1), the cumulative effect of the v/c term, summed 

over 10^° co-rotating stars, becomes substantial. Hence, on RGT, increasing 
rotation velocities are expected in a spiral galaxy. While dark matter might 
exist, there is no need in RGT to hypothesize dark matter just to explain 
non-Newtonian rotation velocities in spiral galaxies. Nor is it necessary to 
make the further hypothesis that the dark matter is peculiarly distributed in 
the form of a halo, just to explain the flattening of the rotation curves. It is to 
be expected on RGT that, at the edge of the galaxy, where the effects of the 
central mass diminish, the velocity-dependent forces between neighbouring 
stars, will make the relative velocities go to zero (i.e., the rotation curves will 
flatten out). Since RGT is based on Lorentz covariance alone, this explains 
the non-Newtonian behaviour of galactic rotation curves in a way which is 
simpler (i.e., makes no hypothesis), and hence to be preferred over both dark 
matter and MOND. 

Since RGT makes no hypothesis, the refutable consequences of RGT are 
not conhned to astrophysics: RGT can be tested in the laboratory. In a 
modihed Gavendish experiment, if the larger (attracting) masses are rotating 
this would result in a torque slightly different from the one when they are 
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static. This also means that static and dynamic ways of measuring the 
Newtonian gravitational constant [13] G may lead to different values. The 
effects of RGT could also be tested by carefully-designed near-earth flyby 
of spacecraft which would be affected by the rotation of the earth as a - 
effect [HI along the lines calculated earlier. 

We emphasize that RGT (which requires only Lorentz covariance) remains 
useful regardless of how it compares with GRT (which demands general co- 
variance). Thus, GRT is simply not usable in practice in situations (like 
the galaxy) where what is used instead is Newtonian gravitation, and RGT 
(since it is Lorentz covariant) is a dehnite improvement on that. Therefore, 
RGT remains the theory of choice for numerous-body problems. Further, 
RGT, rather than Newtonian gravitation, remains the correct flat spacetime 
limit to use. 

1.3 Arguments against retardation 

Nevertheless, there have been objections to retardation in gravity, since pre- 
relativistic times, when naive theories of retarded gravity were hrst proposed 
(to explain the discrepancy in the perihelion advance of Mercury). A key 
objection to those naive theories rests on an argument mentioned by Ed¬ 
dington (which we will call “Eddington’s” argument). The argument is that 
two-body orbits would not be stable with retarded forces. na 

The argument assumes that the two bodies were initially rotating rigidly 
about a common centre of mass. It further assumes that the retarded gravi¬ 
tational force must point to the retarded position of the “attracting” body, 
so that its line of action would not pass through the instantaneous centre of 
mass. On Newtonian mechanics, this means that angular momentum would 
not be conserved. Hence the conclusion that 2-body orbits would be unstable 
unless the “speed of propagation of gravity” is very large, as first asserted by 
Laplace. [T6| This argument does apply to naive theories of retarded gravity, 
as has been argued in the past, im and as we show more rigorously in this 
article. 

However, similar arguments have also been raised against both GRT|18j 
and RGT (in personal communications), and this article aims to show that 
“Eddington’s” argument does not apply to RGT. Now, it is problematic to 
dehne the centre of mass in GRT (or even in special relativity[Tn|). Further, 
there is no physical way for one particle to know the instantaneous position 
of the other, in any relativistic theory. Does that mean that requiring even 
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Lorentz covariance makes all 2-body orbits unstable? In defence of GRT and 
PPN it has been argued, m using Kinnerseley’s “photon rocket”, that the 
velocity-dependent terms in the P^^ cancel the deviation due to retardation, 
so that the net acceleration points to the “instantaneous” centre of mass up 
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to terms of order and the effect of the inexact cancellation shows up only 
in the slowly decaying orbits of binary pulsars. 

1.4 FDEsofRGT 

One can use a similar approach in RGT. However, to proceed in a theoreti¬ 
cally more rigorous way, we first need to take into account a key new feature of 
RGT. The RGT equations of motion Q are retarded functional differential 
equations (FDEs), or state-dependent delay differential equations. Hence, 
the criterion of stability used in “Eddington’s argument” does not extend 
automatically to RGT, since retarded FDEs are fundamentally different from 
ODEs. [12] Even when the retardations are tiny, and seem physically negligi¬ 
ble, it is incorrect to approximate retarded FDEs by ODEs, using a “Taylor” 
expansion in powers of the delay, and this can result in artificial instability, 
no matter how tiny the delay, so long as it is non-zero 

Thus, to examine the problem of stability in RGT in a theoretically rig¬ 
orous way, we first need to solve the FDEs of RGT. 

1.5 Aim 

Accordingly, here we report on a solution of the FDEs of the 2-body prob¬ 
lem in RGT. (The solutions reported earlier for both the galaxy and NASA 
spacecraft were essentially 1-body solutions for a body in the presence of a 
relatively large rotating mass, since the aim there was only to bring out the 
effects of the velocity dependent term in the RGT force, in those situations, 
where the minute movement of the larger body (earth, galaxy) due to the 
smaller one (spacecraft, star) is of no concern and can be neglected.) Our 
aim, in this article, is limited to the question of stability of 2-body orbits in 
RGT, when the full FDEs are solved. This also clarifies exactly how to solve 
the full FDEs of RGT. Other issues will be examined in subsequent articles. 

^For a worked-out counter-example, see m- 
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2 FDEs of motion in RGT 


Since the two equations Q have to be solved simultaneously, it is convenient 
to use a common time parameter, which we take to be the coordinate time t. 
We assume that the functions t = ti(si) and t = ^ 2 ( 52 ) are suitably invert¬ 
ible and (at least) twice continuously differentiable, and will not explicitly 
indicate them further. Thus, we have ^ = 7i, and ;^ = 72 , where 71 and 
72 are the respective Lorentz factors. Using an overdot to denote derivatives 
with respect to t, we have, by the chain rule, ^ = 7iW- 

Similarly, = Tip. 

Hence, ([H can be rewritten 


w 

Ui 


-^ 1 , 

7i 

1 Fi2 
7i rni' 


( 6 ) 


with similar equations for particle 2. 

Since the zeroth component of these equations is not independent, we 
can write them in 3-vector notation using W = (ct, yi(t)), Y 2 = (ct, y 2 (t)), so 
that Yi = (c,Vi), Y 2 = (c,V 2 ). Let ui and U 2 denote the space components 
of the velocity 4-vectors Vi, and V 2 , so that ui = U 2 = 72 'i’ 2 - Further, 
we let r 2 ret = 2 / 2 (^ 12 ) — yiit), denote the 3-vector corresponding to i? 2 ret- 
Then the final equations are 


1 

2/1 = —ui 

7i 

Ul = - fl2 

mi7i 


(7) 


where 


/12 — fl^ 2 ret + bU 2 ret 


( 8 ) 


and 


a = 


kc^ 


V{R2-V2f 


b = ab b 


2 ret 


{R 2 ■ Ul) 


(U 2 • Ul) 


2 ret 


(9) 


The equations ([^ (accompanied by (|^, ([^), together with the corresponding 
equations for particle 2 are the four 3-vector equations (or 12 equations in 
all) we will actually solve. 
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2.1 Past data 


Past data must be prescribed for the solution of the FDEs of RGT (unlike 
the “initial” data appropriate for Newtonian ODEs). We follow our earlier 
strategy of taking the past data to be given by the theoretical 2-body orbits 
of Newtonian gravitation, around the (Newtonian) instantaneous centre of 
mass. Of course, past data could be prescribed in an inhnity of other ways, 
but this strategy suits our immediate purpose, since it allows us to compare 
RGT with Newtonian gravitation (which, we know, gives stable orbits in the 
2 -body case). 

Since our aim is limited to the question of stability, we take up only the 
case of the (mean) Sun-Jupiter system (ignoring other planets). The cor¬ 
responding theoretical, 2-body circular motion is determined in Newtonian 
gravitation by prescribing the mass ratio of Jupiter (=0.000954265748 so¬ 
lar masses), and its mean distance (=5.2 AU). With our strategy of using 
the Newtonian motion as the past data, the same parameters also suffice to 
obtain a solution of the 2-body problem of RGT. 

2.2 Units 

As explained in earlier publications, efficient computational units must be 
appropriate to the past data, and we will use as units, 1 solar mass, 1 AU, 
and 1 earth year. 

2.3 Method of solution 

The technique of solving FDEs has been discussed exhaustively in previous 
publications, |22l |23] and there is nothing new to add. 

The planetary 2-body problem in RGT is not stiff, so the code retard m 
is appropriate. When proceeding on the method of steps[25] (used to solve 
retarded FDEs) the code must be constrained to steps of the size of the delay. 
The delay (~ 8.3 x 10“^ years) is small compared to the time period (~ 100 
years) over which we need to solve the system of equations, to be able to 
resolve the question of stability. Nevertheless, while codes like radar for stiff 
equations can take step sizes larger than the delay, they would be needed only 
in extreme relativistic situations. For the present, using the simpler and long- 
established RETARD code helps to keep the whole calculation transparent. 


3 Results and discussion 


We solved the equations for about 3 centuries, or about 25 Jovian cycles. 
The motion remains stable and the distance between the two bodies differ 
by at most 2 x 10“® AU (~3 km) from the instantaneous distance between 
them in the theoretical (circular) Newtonian motion. The marginal difference 
arises because the orbits do not remain circular but turn into ellipses in 
RGT (as would happen due to any perturbation on the Newtonian theory). 
Had we prescribed elliptic rather than circular past orbits, there would no 
longer be any visible difference between the orbits with RGT and Newtonian 
gravitation. Thus, unlike the galactic case, for the planetary case, RGT agree 
closely with Newtonian gravitation. 

Specihcally, the calculation shows that the delay torque does not induce 
any instability in RGT. We can try to understand this in Newtonian terms. 
The explicit velocity dependent ()() term in the RGT force is essential for its 
Lorentz covariance. On Newtonian intuition, this term makes the net force 
point closer to the “instantaneous centre of mass” and ensures stability. Sup¬ 
pose that term were absent, as happens in “naive retarded gravity” (NRG), 
where the force varies inversely as the square of the retarded distance, and 
points exactly to the retarded position of the other body. In that case of 
NRG, or a “retarded inverse square law”, there is indeed a small instability 

(Fig.§. 

The Lorentz gamma factors are an additional detail to be taken into ac¬ 
count. They were absent in “Eddington’s” argument. These Lorentz factors 
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slow down the velocity of each body by a tiny factor of ^—but, then, in RGT 

(unlike NRG) the delay torque too is at most of that order, for the leading 

order - term cancels. As is well-known, the perihelion shift of Mercury is 
2 

an effect of order but those tiny differences for the case of elliptic orbits 
require a detailed analysis, and we report on them separately. 

Establishing the stability of planetary motion in RGT sets to rest the 
doubt that the higher-than-Newtonian velocities obtained using RGT for 
the case of galaxy were due to instability. Due to the different conhguration 
(of masses), and prescribed past data in the two cases, the explicit velocity 
dependent ()() term, in the RGT force, systematically adds up in the galactic 
case, while, in the planetary case, it cancels to ensure stability. Hence, RGT 
agrees with Newtonian gravity for planetary motion despite giving signih- 
cantly different results for the galaxy. 
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Establishing that the different results for the galaxy are not due to any 
intrinsic instability in RGT opens the way to explore the possibility that 
the accelerating expansion of the cosmos may be similarly due to a different 
conhguration (such as cosmic rotation) rather than yet another invisible new 
substance. 


4 Conclusions 

We have solved the full FDEs of RGT for a model 2-body planetary system. 
The system remains stable, and RGT gives results almost identical to New¬ 
tonian gravity for this case of planetary orbits. Thus, the different results 
that RGT gives for the galaxy are not due to any instability. 
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Comparison of Sun-Jupiter distance 



Time in years 


Figure 1: The relative distance of the two bodies (— ?/ 2 (^))^) stays 
stable. The RGT distance differs by at most 2 x 10“® AU or ~ 3 km from the 
Newtonian distance. This difference arises since RGT changes the prescribed 
Newtonian circular orbit to an ellipse. 
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Comparison of Sun-Jupiter relative speed 



Figure 2: The relative speed of the two bodies (\/(iti — 1 ^ 2 )^) remains stable. 
The maximum difference from the Newtonian case is 10“® AU/year or around 
4.74 X 10“® km/s. 
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Instability of Naive Retarded Gravity 



Time in years 


Figure 3: In naive retarded gravity (NRG) the force varies inversely as the 
square of the retarded distance and points in the direction of the retarded 
position of the other body. The instability of the naive theory is brought out 
by the plot which shows how the relative distance between the two bodies 
iV{yi{t) — y 2 {t)Y) systematically increases over time in the naive theory. 
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